Representing multiqubit unitary evolutions: spin coherences and infinitesimal 

coherences 
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For the tensor of coherences parametrization of a multiqubit density operator, we provide an 
explicit formulation of the corresponding unitary dynamics at infinitesimal level. The main advan- 
tage of this formalism (clearly reminiscent of the idea of "coherences" and "coupling Hamiltonians" 
of spin systems) is that the pattern of correlation between qubits and the pattern of infinitesimal 
correlation are highlighted simultaneously and can be used constructively for qubit manipulation. 
For example, it allows to compute explicitly a Rodrigues' formula for the one-parameter orbits of 
nonlocal Hamiltonians. The result is easily generalizable to orbits of Cartan subalgebras and allows 
to write the Cartan decomposition of unitary propagators as a linear action. 

PACS numbers: 03.65.Ud, 03.67.Mn, 03.67.-a 



I. INTRODUCTION 



The easiest and most promising type of "quantum net- 
work" i.e., of collection of quantum systems to be manip- 
ulated individually or jointly for the purposes of quantum 
information processing, is by far composed of qubits i.e., 
of collections of two level systems. For such systems, in 
we investigated the use of a particular parametriza- 
tion, called the tensor of coherences, obtained by the jux- 
taposition of an affine Bloch vector for each qubit, and 
of widespread use (with minor variations) under differ- 
ent names like cluster operators 0, Stokes tensor Q or 
product of operator basis 0, @ in the literature on NMR 
spectroscopy. It is also closely related to multiparticle 
spacetime algebra and to the recent work of Havel Q • 
Our tensor could be considered a particular parametriza- 
tion of the "nonsymmetric real density matrix" of es- 
pecially suited to emphasize the Lie algebraic point of 
view of the equations of motion. 

The scope of the present paper is to discuss how the 
differential equations describing unitary dynamics must 
be formulated in the tensor of coherences basis. The idea 
that the unitary evolution of a qubit density matrix (pure 
or mixed) given by the Liouville-von Neumann equation 
becomes a linear vectorial ODE for the Bloch vector is 
generalized to multiqubit densities. Mathematically, this 
could be thought of as "passing to the adjoint representa- 
tion" , its starting point being a formula for the decompo- 
sition of nonlocal commutators in terms of local commu- 
tators and anticommutators (see Appendix); practically 
it corresponds still to replacing a conjugation action on 
matrices with a linear action on the vector obtained by 
stacking the columns of the tensor. In particular, when 
operations are local, a unitary transformation reduces to 
a multilinear action, i.e., a linear action on each piece of 
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the tensor of coherences. When instead nonlocal trans- 
formations are used, also their infinitesimal generators 
can be nonsplittable and multispin correlations are in- 
duced. In this case the notation highlights which qubits 
are involved in correspondence of each nonlocal gate. As 
a matter of fact, the major advantage of the formalism 
is that both the pattern of correlations (or "coherences" 
as they are called in the NMR spectroscopy literature) of 
the density tensor and the pattern of the couplings at in- 
finitesimal level become very transparent as both are de- 
composed with respect to the same basis of observables. 
In particular, they both show the same hierarchy of cor- 
relations (that originate from the affine structure of the 
tensors and of the corresponding Lie algebras of genera- 
tors) which allows to keep track of all reduced dynamics 
and reduced densities in a natural way. The idea of as- 
sociating coherences to the degrees of freedom of qubits 
and of manipulating qubits through the corresponding 
Hamiltonians is common for example in the literature on 
spin systems in magnetic fields 0, IH IU E3> El • However, 
the principles apply to any network of qubits. The price 
to pay is a larger dimension of the matrices representing 
the infinitesimal generators: while the size of the Hamil- 
tonians grows as 2™ in the number n of qubits, in the 
adjoint representation it grows as 4™ = 2 2n . 

As an example of the insight gained into the dynam- 
ics of the system, we compute explicitly the integral flow 
of any nonlocal (constant) Hamiltonian, by means of a 
Rodrigues' formula which shows that the exponen- 
tial can be written as a sum of tensor products. Since a 
Cartan subalgebra |13| contains only commuting vector 
fields, also the multiparameter orbit of a set of generators 
belonging to a Cartan subalgebra admits an explicit inte- 
gration. The Cartan decomposition becomes then a con- 
catenation of local and nonlocal linear actions that can 
be expressed directly in terms of the infinitesimal gener- 
ators, rather than of exponentials. Such a decomposition 
has recently attracted considerable attention as a tool for 
constructing universal quantum gates which are optimal 
in the sense of time minimizers or complexity minimizers 
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mm. 

A couple of other examples is discussed, mainly focused 
on the manipulation of qubits in presence of entangle- 
ment. For example we show how to create entanglement 
at distance between qubits that are not directly coupled 
according to two different schemes, one in which the en- 
tanglement is distributed via an entangled ancilla, the 
other via a (always) separable ancilla as in |l5j . 



II. LIE BRACKETS AND ADJOINT 
REPRESENTATION FOR SPIN \ SYSTEMS 

A. One-spin 

Consider the rescaled Pauli matrices and identity ma- 
trix: 



A 



1 



1 o 

1 



Ai 



-i 

1 



1 



1 

1 



1 

-1 



with the commutation relations 

[A , A fc ] = 0, [X u X 2 ] = V2i\ 3 , 
[A 2 , A 3 ] = V2iX u [A3, Ai] = V2iX 2 

and the anticommutators 

{Xj, Afe} = \/25jkXo, 

{Xj, X } = {X Q , Xj} = V2Xj, 



(1) 



j. k G {1, 2, 3}. The operator "ad" is defined as fol- 
lows: adAjAfc = [Xj, Xk] = J2i=o c 'jk^i where opera- 
tions involving the index only produce a null result: 
= 0. Using the "structure constants" 
we obtain an "adjoint basis" associated to the Xj ma- 
trices, given by the four 4x4 matrices ad,\ , ■ ■ • , ad^ 3 of 
purely imaginary entries (ad a 



OA' 



u j0 



■i I kl 



adxn = 



4x4, 



ad Al = v2i 







-1 

.0 1 . 



ad A 



V2i 



0" 
1 

0-100 



ad\„ = \pli 



'0 0- 

0-10 

10 





The Pauli matrices are such that 
form a basis of su(2), while the 
form a basis of so(3) = ad 6U ( 2 ), 



-iX\, — iX 2 and — iA 3 
-iadxj, j = 1, 2, 3, 
the adjoint represen- 



tation of su(2). Concerning the "antiadjoint" operators 
aad^ , j — 0, 1, 2, 3, they can also be defined in the 



same fashion as the adA^, i.e., by means of 4 x 4 ma- 
trices obtained from aad^ Afc = {Xj, Xk} — X^=o s jfeAi 5 
j, k, I <G {0, 1, 2, 3}, so that a linear representation of 
a,a.d\ j is given by (aadAj) fe; = Sj k with the 4x4 matrices 
aadAj easily computed from flfr: 



aad\ 



V2 



a,ad\ 2 = \pl 
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aad a, = \/2 
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B. Two-spin 

Call Aj k = Xj <X> Afc, j, k G {0, 1, 2, 3}. Up to a con- 
stant, the Ajk form the so-called product operator basis, 
see |4j, and are subdivided into spin operators (A o), 
1 spin operators (A i, A 02 , A 03 , Ai , A 20 , A 30 ) and 2 
spin operators (An, A i2 , A i3 , A 2i , A 22 , A 23 , A31, A 32 , 
A 33 ). The set of —iAjk j, k e {0, 1, 2, 3} contains a 
basis of the 9-dimensional tensor product Lie algebra 
su(2) (g) su(2) plus a basis of the 6-dimensional "tensor 
sum" Lie algebras su(2) su(2) arising from the affine 
elements. Just like iXo su(2), so — iA 00 £ su(2) <gisu(2) 
and -iA 00 i su(2) ®su(2). From |jX3j|: 



[Ajk, Aim] = [Aj ® X k , Xi 8) X m ] 

= adA Jfc A/ m = adA j( giA fc A; 

_ 1 

~ 2 
_ 1 

~ 2 



A, 



([Xj, Xi] ® {A fe , A m } + {Xj, Xi} ® [A fe , A m ]) 



ad Aj A; ® aad Afc A m + aad A3 A/ ® adA fc A m ) 



(2) 



In terms of the adjoint representation, J2J can be ex- 
pressed as a 4-tensor, function of the two 2-tensors c l 



and s l jk as: 



ad Ajfe 



adA^igiA 
1 



(ad Aj <X> aadA fc + aad Aj ® ad Afc ) 



' lm 



5(4 



of elements 

( adA i*/(m 2 
so that (j2J becomes: 

[Ajk, Ai m ] — (adAjigiAfc 

-1(4 



km 



^km 



• ' lm A P« 



km 



C lm } A - 



(3) 



(4) 



(5) 
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where we have used the summation convention over re- 
peated indexes (in the range 0-^3). For j ^ and k ^ 0, 
the — iadA Jfc of eq. J3J) form a basis of the adjoint repre- 
sentation of su(2) <g su(2), ad su (2)(g)s U (2) = so (3) (g so (3). 
The remaining elements account for the afhne structure 
i.e., for ad 5U (2)e«iu(2) = so(3) 0so(3). As and s q km are 
4x4 matrices, the resulting Kronecker product adA J( . 
is a 16 x 16 matrix. However, it has a row and a col- 
umn entirely composed of zeros in correspondence of Aoo 
and, given Aj/. with jk ^ 00, $ A; m with (7m) ^ (00) 
such that adAj fc A; m = A o- Furthermore, adA 00 being 
the trivial matrix of all zeros, it is not a basis element 
in the adjoint representation. Also in the adjoint rep- 
resentation the index in a slot corresponds to trivial 
dynamics in the corresponding site. For example 

ad A 3 o = \ (adAj <S> aad Ao + aad A . <g 0) 

1 (6) 
= -^ad A] <g h- 



C. n-spin 

In the n spin case, A.j 1 ...j n — Xj 1 (g . . . (g Xj n , jk £ 
{0, 1, 2, 3}, k £ {1, . . . ,n}, are the basis elements. The 
Lie bracket [A^...^, A^ . fcJ can be computed according 
to the rule IjAlfl . For example, for n = 3 from l|A4(l : 

[Ajki, A mpq ] = i (ad Aj A m <g aad Afc A p (g) aad A; A g 

+aad Aj A m (g ad Afc A p (g aad Al A g 
+aad Aj . A m (g aad Afc X p <g ad Al A g 
+ad Aj . A m (g ad Afe A p (g ad A; X q ) 

= i (ad Aj . <g aad Afc (g aad Ai + aad Aj (g ad Afe (g aad Al 

+aad Aj . <g aad Afe (g ad Al + ad Aj . (g ad Afc <g ad Ai )^* A rst 




(7) 

Remarkably, the building blocks needed for the n-qubit 
case are just the structure constants c l j k and s l j k com- 
puted above. For n spins, the afhne structure prop- 
agates itself throughout and determines a hierarchy of 
subalgebras of tensor product and tensor sum type. 
The — iAjj.j^ , {j\ ■ ■ ■ j n ) 7^ (0...0), form a joint ba- 
sis of the Lie algebras su(2)®", su(2) su(2)^(™" 1 ), 
. . ., su(2) © ... su(2) (plus all factor permutations) 
and -iadA^ . n , {ji . ■ ■ j n ) ^ (0...0), a joint basis of 
ad S u(2)®": a d sl i(2)esii(2)®c™-i)) ad 5U (2) e ... e5U (2) (plus, 
again, all factor permutations). In both notations, the 
number and position of the indexes "0" uniquely deter- 
mine which spins are involved into the — iadA . 



III. UNITARY EVOLUTION IN TERMS OF 
THE TENSOR OF COHERENCES 

For qubits, the same basis elements Aj 1 ...j n that de- 
scribe the infinitesimal generators can be used also for the 
density operators. This is well-known in the literature on 
spin systems under the name of "coherences" , Q , and can 
be formalized in terms of 4 x 4 x . . . x 4 tensors which we 
call tensors of coherences. See 0, 0, 0, El 03 E f° r an 
overview. The scope of this Section is to show how ten- 
sors of coherences and adjoint representation fit together 
in the description of unitary dynamics of multiqubit den- 
sities. 



A. Density operators and tensor of coherences 

This Section follows [l|. The Aj 1 ...j n form a complete 
orthonormal set for Hermitian matrices and can be used 
to obtain an afhne tensorial representation of the density 
operator of n qubits: p = Q ix "^ n Aj l ...j n ,jk £ {0, 1, 2, 3}, 
k E {1, . . . , n}, with g? 1 ■■■■ J ™ = tr {pAj lt „j n ) the expecta- 
tion value along the observables Aj liti j n . This represen- 
tation has several advantages briefly recalled below: 

• it captures all degrees of freedom of a density op- 
erator; 

• each term g-? 1 - ^™ in the tensor depends on a cer- 
tain number of qubits: this is uniquely determined 
by the number of nonzero indexes in the sequence 
ji . . . j n . The pattern of nonzero indexes also iden- 
tifies which qubits are involved, formalizing the idea 
of coherences of widespread use for spin systems. 

• all correlations of all orders and all reduced den- 
sities are already contained in the tensor: tracing 
out a qubit means collapsing the corresponding in- 
dex to and rescaling everything by \/2, For ex- 
ample, if pa 2 ...a„ — t r Ai (p) = g- 72 '"- 7 ™ Aj 2 j n then 

• Since 

tr(A jfc A im ) = Sji6 km , (8) 

j, k, I, m G {0, 1, 2, 3}, the degree of mixing be- 
comes the Euclidean norm of g 31 "'■ J ™: 

tr(p)= ]T (e*-*.) a (9) 

ji—j n =0 

and hence, since q° - ° — (l/y/2) n , for (ji . . . j n ) ^ 
(0...0) t he tensor £ Sf" 1 G R 4 "" 1 with 

< r ^ - (e°-°) 2 = \/l-(l/2f; 

• complete mixing corresponds to r = (i.e., to the 
null tensor except for the afhne constant g°—°): 
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• pure states correspond to r = y/l — (1/2)™; 

• uncorrelation corresponds to gj^—Jn = 

tflA ■ ■ ■ where e?Xi = {V2) n ~ 1 Q jl0 -° is the 
4- vector of the reduced density pa 1 — trx a ...A„ (p) 
and so on; 

• partial transposition of a qubit becomes a change 
of sign in the terms having index 2 in the corre- 
sponding slot. For example 



p T Al = Q Oh..dn Ao , 2 ___. n+e lh..Jn Alh ___, n 



(10) 



and so on; 



• checking bipartite entanglement can be done by the 
simple test (fTT)|) . 



B. Liouville-von Neumann equation 

The Liouville-von Neumann equation for the n-qubits 
density p is: 



P = -i[H, p] = -iadn(p) 



(11) 



where H = H' is the Hamiltonian of the system. From 
Section [HJ we have that H = ft,- 71 '" 3 " A^ j , jk £ 
{0, 1, 2, 3}, k£{l,...,n}. 

If we have two qubits then, in terms of the tensor of 
coherences, eq. <|11[) corresponds to: 



ih jk 



(<5 



, <8>s? +s p ,®cl )p lm . 

I km 1 il km I = 



(12) 



In order to show (|12f> . derive g P9 = tr (pA pq ) and use (QJ 
and JEJl: 



= tr(pA M )=trH[ff,p]A M ) 
= tr(-i^[A^, A /m ]f/ m A P9 ) 



( Cj r , 4 m + sj, ® c s km ) g lm tr (A rs A pq ) 



The component of the Hamiltonian along Aoo is irrele- 
vant: even if h Q0 ^ it must be -ih 00 &d Aoo = 0. The 
meaning is similar to the single spin case: global phases 
are neglected in (|11[> and l|12l) . 

Since 112(1 is a linear system, if h 3 are constant the 
integration can be carried out explicitly: 

gP1{t) = ( e -* thJk ^ ]k Y q g^(0). (13) 



Notice that when 2-spin coherences are lacking, h? = 
V j ^ and k ^ 0, i.e., when only LOCC operations 
are performed, the exponential in (|13fl splits. In fact, 
[Ajo, Aok] — and therefore the infinitesimal generators 
Ajo and Aofc can be "reduced" as well. From ©, the 
unitary propagator in eq. i|13[l becomes: 



\/2 A j 



0* 



it(W°ad A . +h ofe ad Aofe ) = /-«fc*>»d Aj0 ^ 
h) (h 



V5 



where the factor A= comes from (0. Therefore 



e _ it iL_ ad> 



1 

S*0(3) 



1 

SO(3) 
(14)' 

which allows the state to evolve on at most a 6-parameter 
orbit sitting inside the 15-dimensional affine sphere §* 5 , 
with r identified by JSJ. If p(0) is separable then under 
(|14fl so is pit) for all t, the 6-dimensional manifold con- 
tains all the separable states. When instead the Hamil- 
tonian has h? 7^ for j ^ and k ^ 0, the evolution of 
the two qubits becomes coupled. 

Similarly to the 2-qubit case, if we have n qubits we 
obtain: 

? p >ii j j ( ad An..,Jt:t ^'" fe "' 

where adA- is computed as in Section III CI 



IV. INTEGRAL FLOW OF NONLOCAL 
HAMILTONIANS 

We first restrict to 2 qubits, although all arguments 
generalize to n qubits. First we give an explicit for- 
mula for the integral of each "elementary" generator 
Ajk ■ From Section III Al we have that ad A . aad A . = 
aad Ai ad Ai = 0. This implies that the series expansion 

exp (-iiad Ajfe ) = J2p°=o ^T~ ad A jfc nas a particularly 
simple expression, since for all p 

ad a = — fadx ® aad, + aad^ <g> &<S 

Aj-Jj 2P \ A 3 A fc A 3 A ' 

The powers of ad Aj and a,a.d\ j are easily computed since 
ad Aj and aad Aj are diagonal and "complementary" : 

• if j = 1, ad Ai = 2(5 3 3 + fo), aad Ai = 2(6 U + S 22 ); 

• if j = 2, ad A2 = 2(5 2 2 + M, aad^ = 2(S U + S 33 ); 

• if j = 3, ad Aa = 2(5 2 2 + $33), aad Aa = 2(6 u + 644); 

so that ad Aj + aad Aj = 2I4. Cubic powers instead are 
ad A . = 2ad A . and aad A . = 2aad A , hence ad A . = ad A k - 
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We can therefore explicitly write down the sum of the 
series as 



exp (— it&dA jk ) = I4 ® I4 — i I t 



2! + 4! 



i 3 t* 

3! + 5! 



ad a 



ad A , 



or, adding and subtracting ad Ajfc , 

exp (— itadA j.) = h<S>h — isin(t)ad\, k — (1 — cos(i))ad A . fc 

(15) 

where the extra terms added are needed because the zero 
order terms do not match: / 4 ® / 4 ^ ad 2 ^ ■ Notice that 
a formula equivalent to l|15(l was used for the same pur- 
poses as ours in Q. Both are tensored versions of the 
Rodrigues' formula for rotations, see 0], p. 291 or [20| . 
p. 28. The splitting is into skew-symmetric (— ia,d\ k ) 
and symmetric part (I4 (g) I 4 and ad^. fc ) of the flow [24j| . 
Notice that both — iadA jfc and ad 2 . fc are tensor products 
of matrices. The nonlocality of the Hamiltonian of i|15|) 
reflects in the fact that we do not obtain a "single" tensor 
product but rather a sum |2j| . Clearly the overall evolu- 
tion of (|15|) is orthogonal. However, the single pieces do 
not describe rotations, neither locally nor globally. 

The same argument can be repeated for any num- 
ber of qubits. For example for 3 qubits we have 

exp (-itadA jkl ) = E^Lo iz ?- ad A Jfc! > with 

a,d A , k = i (ad^. <8 aad^ fc (g) aad^-|- aad^ ad^ fc (g) aad^ 

+aad^ (g) a ad^ fc (g) ad^ ; + ad^. (g) ad^ fc ® ad^ ^ 

'(16) 

where now ad a = A ad a ., , . The sum of the series is 
then 



exp (-ita,d Ajkl ) =lf 3 - iv / 2sin(-^=)ad Ajfci 



(17) 



So far we have only considered a single "coordi- 
nate direction" (Ajk for the 2-qubit case). The for- 
mulae however extend in a straightforward manner to 
linear combinations of commuting generators, even de- 
pending on more than one parameter. A maximal 
orbit of integrable flow is obtained obviously in cor- 
respondence of a Cartan subalgebra [l3l Il4j . i.e., a 
maximal commuting subalgebra in the Lie algebra of 
nonlocal operations of the system. For the 2-qubit 
case, let us concentrate on the "nonlocal subalgebra" 
ad BU ( 2 ) ( g ); i U (2) = so (3) <8>so(3). A Cartan subalgebra is for 
example given by f) = span{— iad All , — *adA 22 , — iadA 33 } 
(or by span{— iadA 12 , — iad Aai , — iadA 33 }, etc.). The 3- 
parameter orbits of such subalgebras are integrable as 



can be seen by the splitting of the exponential 
exp (—i (a 11 adA 11 + a 22 adA 22 + a 33 adA 33 )) = 

= exp (— ia 11 adA 11 ) exp (— ja 22 adA 22 ) exp (— ia 33 adA 33 ) 

(18) 

for real o?K The "marginal" subalgebra of local oper- 
ations 00(3) © 00(3) does not commute with the Cartan 
subalgebras. It is known that [00(3)000(3),!)] gener- 
ates the entire 15-dimensional Lie algebra 00(3) 000(3) U 
so(3)(80o(3) and that "exponentiating" this splitting we 
get the Cartan decomposition of the corresponding Lie 
group. With our formalism, such conjugation action be- 
comes a linear action, obtained by the concatenation of 
a bilocal exponential as l|14|) and of (|18fl . In other words, 
any unitary operation acting on a 2-qubit density can be 
written for the tensor of coherences as the concatenation: 



-ad;. 



v'2 



ad,\ 



exp (- 



,11 



ad All/ 



exp (— ia 22 adA 22 ) exp (— ?o; 33 adA 



for real a? k . Each exponential can be replaced by the cor- 
responding sum of tensors (given by (|15fl for the nonlocal 
pieces and by exp (— it&d\ j ) = I a — -4^ sin(v / 2i)ad> j — 
(i-cos(V2t)) ^ ^ e one-parameter orbit of a single 
qubit). 



V. EXAMPLES 

In Example IV Al it is shown how to express in term of 
the tensor of coherences the discrete unitary propagator 
corresponding to a standard 2-qubit gate, the C-NOT 
gate. In the three qubits of Example IV Bl entangling be- 
tween two "distant" qubits is achieved through indirect 
coupling by means of an entangled ancilla. In Exam- 
ple IV CI instead, for the same purposes the scheme of 
15] is used, in which the ancilla remains separable for all 
times. 



A. C-NOT gate 

It is well-known that elementary gaits of a quan- 
tum computer being discrete unitary operations, they 
can written in terms of the corresponding infinitesimal 
Hamiltonians. In particular, in the literature on quantum 
information processing by means of NMR spectroscopy 
|5( this was done in terms of the product of operator 
bases of which our formalism is just a variation. For ex- 
ample, in correspondence of the computational basis of 
two qubits 1 00), |01), 1 10) , 1 1 1) , the Hamiltonian of the 
C-NOT gate 



U c - 



NOT — 



1 0' 

1 

10 

10 



6 



is given 



by i 



B. Three-qubit: entangling at distance (I) 



H, 



C-NOT 



o o ■ 

10-1 


0-101 



In terms of the Ajk, this is -Hc-not 
\ (A o - A 3 - Aio + Ai 3 ), and therefore for g jk 
have the orthogonal matrix 

-Rc-not = e~'^( adAno ~ adAo3 ~ adAlo+adAl 3), 
which computed by means of JSJ yields 



wo 



-Rc-NOT 



- 1 
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|00) 
|01) 
|01) 

111) 



If we are given the 4 computational basis states 

1/2,0,0,1/2,0,0,0,0, 
0,0,0,0,1/2,0,0,1/2 

1/2,0,0,-1/2,0,0,0,0, 
0,0,0,0,1/2,0,0,-1/2 

1/2,0,0,1/2,0,0,0,0, 
0,0,0,0,-1/2,0,0,-1/2 

1/2,0,0,-1/2,0,0,0,0, 
0,0,0,0,-1/2,0,0,1/2 



it is straightforward to check that i?c-NOT behaves as 
a C-NOT gate with the second qubit acting as control 
qubit. Notice that i?c-NOT is not traceless, hence we 
have an Hamiltonian with h 00 ^ 0. As mentioned above, 
this is irrelevant because adA 00 = 0, i.e., in the ad- 
joint representation one always obtains the corresponding 
traceless Hamiltonian. 

The structure of the basis used indicates that -ffc-NOT 
is a non-local operation since it contains A13 (and the 
splitting into basis elements is obviously unique). While 
it leaves unentangled the computational basis elements, 
the same is not true in general for any state. 

Comparing Uq-not and i?c-NOT, the price to pay 
in order to use the tensor of coherences parametrization 
is a larger dimension of the operator involved. On the 
other hand, the matrices are normally sparse and the 
formalism allows to perform the same operation also on 
mixed states. 



Assume we have available coupling Hamiltonians be- 
tween A and B and between B and C. B can be thought 
of as an ancilla being first entangled with A and then 
sent to interact with C. Given a state in which A is max- 
imally entangled with B while C is separable from the 
two (and known), we want to transfer the entanglement 
from the pair (AB) to the pair (AC) leaving B unentan- 
gled at the end of the evolution, without making use of 
coupling Hamiltonians between A and C. Assume pab(0) 
is in the pure maximally entangled state 



{00,11,23,32} 



Q 



(0) 







otherwise. 



and PC = + ^i)- The desired task is accomplished 

in half of the period t p — for example by the fol- 

lowing piecewise constant Hamiltonian: 



-iaAiiit) = 



-iadA 033 
-iadA 220 



t G [0, T -f) 

i (C 112. l£ 

^14' 2 



We obtain also that Pab(0) — PAci^/V) and pb(tt/2) = 
pc(0), see Fig^ As can be seen from Fig. El at the 
entanglement swaps from the pair AB to the pair AC. 
The scheme can be iterated to n qubits. 




FIG. 1: The 64 components of p versus t. 



C. Three-qubit: entangling at distance (II) 

While the previous example is rather straightforward, 
in the literarture there exist more sophisticated and 
surprising methods to distribute entanglement. In [l5| 
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eigenvalues of p and o( p T A 



igenvalues o( p and o( p T E 




2 3 4 



eigenvalues o( p and o( p T c 




FIG. 2: The eigenvalues of p (dashed lines) and of the 3 partial 
transposes p TA , p TB and p T ° (solid lines) versus t. 



it is shown that for the 3-qubit separable state p m — 

\ (ELo \M-k, o) (i> k , i)- k ,o\ + E)= li. 3, i) (i, j, i 



with |^ fc ) = (|0) + e tk */ 2 \l}) /y/2, it is possible to find a 
cascade of two C-NOT gates, one with C as control qubit 
and acting on A and the other with B as control qubit 
and acting on C, such that at the end of the operation 
A and C are both entangled but for the whole process 
B remains unentangled. In terms of the Hamiltonian of 
the C-NOT computed in Section IY Al this is equivalent 
to the following piece-wise constant 3-qubit infinitesimal 
generators, obtained permuting the indexes of -ffc-NOT 
and adding a "0" in the correct slot [2|J: 



-iad#(t) 



—% (— adA 300 
—i (— adA 003 



ad Aool +ad A301 ), t€[0,-J=) 
adA 010 + adA 013 ) , t £ [-^=, 



If X = 672' th6n 



1 



Pin 



Pint = 



2^2 
-2A220 
1 



Pfin 



Aooo + 2A003 + 2A110 + xA n3 
xA 22 3 + xA 330 - xA 333 , 
Aooo - 2A033 + xAin - 2A122 
xA 22 i + xA 30 3 + xA 330 , 
xAioi + xAi3i 
xA 23 2 + 2A303 + xA 333 , 



272 

-xA 2 12 

-— ^=A oo — a:A 3o 
272 

-a;A 20 2 



where p- ln t is the density after the first C-NOT gate and 
pan the final state. Simulating the evolution of the sys- 
tem, we get that indeed B maintains a positive partial 
transpose (PPT) for the whole interval, as can be seen 
in Fig. |3 while A acquires a negative partial transpose 
(NPT) in the first half and keeps its through the second 
half. In this second part also C shows NPT. The behav- 
ior can be explained in terms of bipartite entanglement 
of different cuts of the 3 qubits. Compare Fig. |3| and 



Fig. ^] 27]. In the first half of the interval, A is entan- 
gling itself with the 2-qubit reduced density psc- Such 
entanglement is bipartite and is not "visible" at the level 
of 1-qubit reduced densities of B and C. The same thing 
happens between C and (AB) in the second half of the 
operation. The example is a well-cooked one as for all 
times there is no entanglement showing between B and 
(AC) (not just "at the end" of the gate). The doubt 
that remains is whether the final result is truly creation 
of entanlgement between A and C, or rather is only a 
"superposition" of two 1-qubit - 2-qubit bipartite entan- 
glement. Notice that a third C-NOT operation on A 



eigenvalues of p and of p a 



eigenvalues of p and of p T B 



eigenvalues of p and of p T c 



2 3 4 



FIG. 3: The eigenvalues of p (dashed lines) and of the 3 partial 
transposes jO Tj4 , p TB and p T ° (solid lines) versus t. 



eigenvalues of p and of p a 



eigenvalues of p and of p 



eigenvalues of p and of p 



2 3 4 



FIG. 4: The eigenvalues of p (dashed lines) and of the 3 partial 
transposes p Tsc , p TAC and p TAB (solid lines) versus t. 



and C (with either of the two as control qubit) leaves all 
three qubits with PPT. 
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APPENDIX A: FORMULAE FOR LIE BRACKETS 
OF TENSOR PRODUCT MATRICES 

Proposition 1 Given Ax, . . . , A„, Bx, . . . , B n E M m , 

the commutator of Ax ® ... ® A n and Bi ® . . . ® B n is 
given by 



[Ax ® . . . ® A n , Bx ® ■ ■ ■ ® B n ] = 
= Ei((4 Bx)®(A 2 , B 2 ). 



(A n , B n )) 

(Al) 



where in each summand the bracket ( • , ■ ) is 

{[•,•] k times, k odd 
{ • , • } n — k times 

and the sum is over all possible (nonrepeated) combina- 
tions of [ • , • ] and { • , • } and over all odd k € [1, n\. 

The anticommutator of A\ ® . . . ® A n and B\ ® . . . ® B n 
is given by 



{Ax ® . . . ® A n , Bx ® . . . ® B n } = 
= E^l(( Al ' B i)®(^2, B 2 )® 

where in each summand the bracket ( 



(A n , B n )) 

(A2) 



is 



{■>■} 



k times, k even 
n — k times 



and the sum is over all possible (nonrepeated) combina- 
tions of [ • , • ] and { • , • } and over all even k £ [1, n\. 

Proof. We will prove the Proposition by induction. The 
formula i|Alfl is obviously true for n = 1 (for n = 2, 3, 4 



it is explicitly given below). Assume it is true for n — 1 
and write a — A x <£> . . . <g> A n _i, f3 = Bx® . . . (S> B n -x- 
Then for n we have 



[a <gs A n , 13 ® B n ] =a(3® A n B n - [3a® B n A n 
± - (a/3 <g) B n A n + (3a® A n B n ) 

= i ([a, (3} ® {An, B n } + {a, [3} ® [A n , B n ]) . 



If [a, [3] contains an odd number of commutators, so does 
[a, 0\® {A n , B n }- Likewise, if {a, (3} has an even num- 
ber of commutators, {a, f3}®[A n , B n ] has to have an odd 
one. If [a, (3] and {a, (3} contain all possible nonrepeated 
combinations of commutators and anticommutators, so 
does the expression [a ® A n , (3 ® B n ], and the induc- 
tion is thus completed. Concerning the anticommutator 
(|A2|) . the same induction arguments can be repeated for 
the following expression: 



{a ® A n , (3 <g> B n } = a(3® A n B n + Pa® B n A r , 
±^(a/3® B n A n + f3a® A n B n ) 



- ([a, 13] ® [A n , B n ] + {a, 0} ® {A n , B n }) ■ 



□ 



While we are not certain of the complete novelty of the 
recursive formulas (|Alj) and (|A2|1 . we are sure that var- 
ious equivalent variants of them are well-known |28| for 
low-dimensional tensors. Restricting to recent related lit- 
erature, check for example [y, \13L l21l | . The commutators 
for the first cases used in the paper are given explicitly 
below. 



[Ax ® A 2 , B Y ® B 2 ] = A Y B X ® A 2 B 2 - B\A\ ® B 2 A 2 

= \ ([Ax, Bx] ® {A 2 , B 2 } + {£?!, Ai} ® [A 2 , B 2 ]) , 



(A3) 



[Ai ®A 2 ® A 3 , B 1 ®B 2 ® B 3 ] = A\B\ ® A 2 B 2 ® A 3 B 3 - B\A\ ® B 2 A 2 ® B 3 A 3 

= 1 ([Ai, B{\ ® {A 2 , B 2 } ® {A 3 , B 3 } + {A l5 Bx} ® [A 2 , B 2 ] ® {A 3 , B 3 } 
+{Ai, Bx} ® {A 2 , B 2 } ® [A3, B 3 ] + [Ax, Bx] ® [A 2 , B 2 ] ® [A 3 , B 3 ]) , 



(A4) 



[Ax ®A 2 ®A 3 ® Ax, Bx® B 2 ® B 3 ® B 4 ] = A X B X ® A 2 B 2 ® A 3 B 3 ® A 4 B 4 - B X A X ® B 2 A 2 ® B 3 A 3 ® B 4 A 4 
= i([Ai, B{\ ® {A 2l B 2 } ® {A 3 , B 3 } ® {A 4 , B 4 } + {A 4 , Bx} ® [A 2 , B 2 ] ® {A 3 , B 3 } ® {A 4 , B 4 } 

+ {Ax, Bi} ® {A 2 , B 2 } ® [A 3 , B 3 ] ® {A 4 , B 4 } + {Ax, Bx} ® {A 2 , B 2 } ® {A 3 , B 3 } ® [A 4 , B 4 ] 
+ [Ax, Si] ® [A 2 , B 2 ] ® [A 3 , B 3 ] ® {A 4 , B 4 } + [Ax, Bx] ® [A 2 , B 2 ] ® {A 3 , B 3 } ® [A 4 , B 4 ] 
+ [Ai, Bx] ® {A 2l B 2 } ® [A 3 , B 3 ] ® [A 4 , B 4 ] + {A 4 , B{\ ® [A 2 , B 2 ] ® [A 3 , B 3 ] ® [A 4 , B 4 }). 



(A5) 
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Even the sign difference with respect to the standard 
SO(3) formula is due to the fact that here the skew- 
symmetric generator is — iadA jfc • 

This does not mean that we have separable superoper- 
ators p2l however, since unitary operators yield "pure" 
quantum operations |23j|. 

Notice that the time interval is rescaled with respect to 
the 2-qubit case of Section IV Al because of the effect of 
the third qubit, see @. 

Obviously Fig. 0] is totally redundant with Fig. [3] 



(P ) =P 
And trivial 
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since it is enough to replace AB 



2 ([A, B] + {A, B}) in the brute force calculation of the 
commutator/anticommutator and regroup appropriately. 



